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Abstract. In his work studying the Zeta functions of famihes of hypersur- 
faces, Dwork came upon a one-parameter family of hypersurfaces (now known 
as the Dwork family). These examples were not only useful to Dwork in his 
study of his deformation theory for computing Zeta functions of families, but 
they have also proven to be extremely useful to physicists working in mirror 
symmetry. A startling result is that these families are very closely linked to hy- 
pergeometric functions. This phenomenon was carefully studied by Dwork and 
Candelas, de la Ossa, and Rodriguez- Villegas in a few special cases. Dwork, 
Candelas, et.al. observed that, for these families, the differential equation as- 
sociated to the Gauss-Manin connection is in fact hypergeometric. We have 
developed a computer algorithm, implemented in Pari-GP, which can check 
this result for more cases by computing the Gauss-Manin connection and the 
parameters of the hypergeometric differential equation. 

1. Introduction 

In his work studying the Zeta functions of famihes of hypersurfaces, Dwork came 
upon a one-parameter family of hypersurfaces in P"~i (now known as the Dwork 
family), defined by: 

Xx : Xi -\ \- x'^- nXxi • • • x„ = 0. 

These examples were not only useful to Dwork in his study of his deformation 
theory for computing Zeta functions of families, but they have also proven to be 
extremely useful to physicists working in mirror symmetry (c.f. [5] ). 

A startling result is that these families are very closely linked to hypergeometric 
functions. This phenomenon was carefully studied by Dwork in the cases where 
n = 3, 4 (see for example [10]) and for n = 5 by Candelas, de la Ossa, and Rodriguez- 
Villegas ([7], i). 

In particular, Dwork's work gives a construction of modules isomorphic to the 
middle (relative) deRham cohomology, equipped with an integrable connection 
which is equivalent to the Gauss-Manin connection. By the work of Katz and 
Oda [19], we know that this connection is essentially differentiation of cohomology 
classes with respect to the parameter. Differentiating each basis element in the 
module gives us a first-order system of differential equations. Dwork, Candelas et. 
al. noticed, in the examples they studied, that the differential equations obtained 
through this method are hypergeometric differential equations. 

In recent years, Dwork's ideas have been generalized to compute Zeta functions 
using p-adic and £-adic cohomology. In studying the Zeta function using £-adic 
cohomology, Katz proved that there was a link between more general monomial 
deformations of Fermat hypersurfaces (of which the Dwork family is an example) 
and hypergeometric sheaves [18]. Rojas-Leon and Wan, independently from Katz, 
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implemented the same approach to compute Zeta- functions ^23^ . In ^20, , Klooster- 
man shows that the p-adic Picard-Fuchs equation associated with the Dwork family 
is hypergeometric. 

In this paper, we follow the direct approach originally used by Dwork and then 
Candelas et al, rather than the rigid cohomology approach. That is, we use Dwork's 
original construction of a module over C and algorithmic methods based on combi- 
natorics and linear algebra. . We have developed computer algorithms, which have 
been implemented in Pari-GP pTj.The GP scripts can be found in the Appendix. 
One algorithm computes the matrix for the Gauss-Manin connection associated 
with the Dwork family, by blocks. Another uses a block of the connection matrix 
to compute the parameters of the associated hypergeometric differential equation 
(and in the process proving that this differential equation is hypergeometric). 
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[24j , she would like to primarily thank her advisor Fernando Rodriguez- Villegas for 
his guidance, support, and great ideas. Many people contributed to the progress 
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Frits Beukers, and Daqing Wan, and the Arizona Winter School for allowing her to 
meet these great mathematicians. Finally, the author would like to thank Michelle 
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and Jonathan Webster for his computational number theory advice. 

2. Background 
2.1. Hypergeometric Functions. 

Definition 2.1. Let, A,B G Z and ai, ... ,aA, Pi, Pb G Q, with all of the 
Pi > 0. The generalized hypergeometric function is defined as the series 
(taking z E C) 

AFB{ai,. . . ,Q!a;Pi, ■ • ■,Pb\z) = 2^ 



(/3i)fc • • • (Mfcfc! ' 



where we use the Pochhanimer notation 



{x)k = x{x + l)---(x + k-l)^ ^-fe^- 

r(a;) 

The Ui will be referred to as "numerator parameters" and the Pi as "denominator 
parameters" . 

Sometimes we will use the shortened notation 

AFB{a]P\z) = AFB{ai, . . ■ ,aA] Pi, ■ ■ ■ , Pb\z). 

Let 9 denote the operator z — . The series aFb^ch] P\z) satisfies the differential 

dz 

equation 

{e{e + pi-i)---{e + Pb-i)- z{e + ai)---{9 + a^)} y = o. 

Following the notation in [2j, 



D{ai, ...,aA;Pi,...,PB)^e[e + Pi-l)---{e + PB-l)- z{e + ai)---{0 + aA). 
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If A = B + 1 this is a Fuchsian differential equation with regular singularities at 
z — 0,1,00 (in Section [5] we will review these definitions). We will focus only on 
hypergeometric functions with this property. 

Notice that the parameters ai , /?; completely characterize the hypergeometric 
function and its corresponding differential equation. 

2.2. Hypergeometric Groups. In Section [5l we will use a certain property of 
monodromy groups in order to relate the Gauss-Manin connection to hypergeomet- 
ric functions. First, we need some definitions from [2 . 

Let H be the fundamental group 7ri(P^ \ {0, l,oo},2o) where zq is some fixed 
base point, for example zq = 5. Then clearly H is generated by go,gi,goo with the 
relation (jooSiffo = 1, as pictured below. 




5oc. 



Figure 1. The generators of 7ri(P^ \ {0, 1, 00}, zq) 

Recall that the differential equation for a hypergeometric function of the form 
nFn^i{a; (3\z) is Fuchsian with regular singular points 0,1, 00. Around a regular 
point, for example zq = i, there are n linearly independent analytic solutions with 
a non-zero radius of convergence. Let A, B,C € GL{V) be determined by analytic 
continuation of solutions along the generators of 7ri(P^ \ {0, 1, 00}, zq), so that 

A ^ goo 
B ^ go 
C ^ gi 

The group F C GL{V) generated by A,B,C with the relation ACB = Id is 
called the monodromy group, and the map 

H GL{V) 
5oo,go,5i A,B,C, 

is a representation of H. 

Definition 2.2. Let y be a finite dimensional complex vector space. A linear map 
g G Gh{V) is called a reflection if 17 — Id has rank one. The determinant of a 
reflection is called the special eigenvalue of g. 
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Definition 2.3. Suppose ai, . . . , a„, &i, . . . , fe„ G C* with aj ^ bk for all j — 
1, . . . ,ri. A hypergeometric group with numerator parameters ai, . . . ,an and de- 
nominator parameters 6i, . . . , 6„ is a subgroup of GL(n, C) generated by elements 
ho, hi, hoc G GL(n, C) such that hochiho = Id, 

n 

det(z - hoc) = ]^(z - Cj) 

i=l 
n 

det{z~h^') = l[{z~bj), 

and hi is a reflection in the sense of Definition 4.1. 

Then we have the following useful result. 

Proposition 2.1. [2] Suppose ai, . . . , a„, 6i, . . . , &„ £ C* wiit/i ^ bk for all 
j,k = 1, . . . ,n and assume bn = 1. Let ai, . . . , a„, . . . , /3„_i G C 6e such that 
Oj — e^^^°'^ for j — I,. . . ,n and bk = e^'"^^*' for k — 1, . . . ,n ~ 1 . Then the 
monodromy group of the hypergeometric equation 

D{ai, ...,«„; ^1, ... , I3n~i)y = 

is a hypergeometric group with parameters ai, . . . , a„, foi, . . . , 6„. 

The most important consequence of this result is that if we have a hypergeometric 
group r in the sense of the Definition 2.3, we can find a hypergeometric differential 
equation whose monodromy group is F. 

3. Dwork's construction 

Suppose X is a hypersurface defined over P", and so it is n — 2 complex dimen- 
sional. The middle deRham cohomology will be the {n — 2)-th cohomology. It is a 
classical result by Lefschetz that the i-th deRham cohomology of a projective hy- 
persurface of dimension n — 2 is identical to the j-th deRham cohomology of P"~^, 
for i 7^ n — 2. In other words, the middle cohomology is the only "interesting" one. 

Recall the Dwork family of hypersurfaces defined by 

: .t" H 1- a;" - nXxi • • • x„ = 0. 

It is not difficult to see that X\ is not smooth if and only if A is an n-th 
root of unity. Let T — C — fi„. It follows that X\ is non singular for A G T. 
Dwork constructed modules over C isomorphic to the relative deRham cohomology 
H^^^{Xx/T), which are quite combinatorial in nature. 

To give Dwork's construction, we will use the notation established in [T^. Let 

F\{xi, . . . ,Xn) := Xi -\ hx" —nXxi • • • a;„. Let £ be the free module (over C(A)) 

generated by the monomials 

. . . T™" — T™ 

with all the Wi > and Wi = mod n. 

Let 2^ be the submodule generated by monomials x'^ with all Wi > 1. Let Di 
be the mapping defined by 

D,: 2^ 2, 

So D^ix"") = w.x"" +nx'lx 



BP 

D.ix"-) = w,x^ + x^^x- 
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We now define the C(A)-module 

Thus, W is a vector bundle over T. 
Notation. We will frequently represent a monomial by its exponent w. 

Algorithm 1 (Reduction Algorithm). We define an algorithm on n-tuples w using 

the relations given by the quotient described above. The input for this algorithm is 
c * w = c * {wi,W2, . . . , Wn), where c is an element of C(A), and w is an n-tuple 
representing a monomial in >V where Wi > n for some i. The output is a list 
that represents how to write cw as a linear combination of monomials for which all 
exponents are less than n, thus "reducing" w. 

(1) Initialize two empty lists, L and M. 

(2) Let i be the first entry such that Wi > n. Define u — cX * (wi + 1,W2 + 
l,...,Wi-n+l,...,Wn + l) and V = c^^^^^ * {wi, . . . ,Wi - n, . . . ,Wn). If 

n I 

(a) If Vi < n for all i, append v to list L, unless {wi, . . . , w;,; — n, .... (i'„) 
is already in the list, then add c "^'"' to the existing coefficient. Go to 
step 3. 

(b) Ifvi > n for some i, append v to list M , unless (?i>i , . . . ,Wi — n, .... Wn) 
is already in the list, then add c "^™' to the existing coefficient. Go to 
step 3. 

(3) If Ui < n for all i, append u to L and proceed to step 4 (here if u is already 
in L we add coefficients as before). Otherwise, proceed to step 5. 

(4) If the list M is empty, output the list L. Otherwise, proceed to step 8. 

(5) If u = w, proceed to step 6. Otherwise, proceed to step 7. 

(6) If the list M is empty, output the list L but with all coefficients divided by 
1 — A" . Otherwise, proceed to step 8. 

(7) Append u to M and proceed to step 8. 

(8) Let a be the first element in M. Apply step 2 to a instead of w and remove 
a from M. 

Proof. Notice that we stop the algorithm whenever the list M is empty. The list 

M contains monomials with entries that arc greater than or equal to n. To prove 
the algorithm terminates we need to prove that using this reduction process we can 
always empty the list. 

Given any starting n-tuple w, the worst that can happen is that every entry 
is greater than or equal to n. In this case, the first step would be to set u = 
cA * {w\ — n + 1, W2 + 1, • • • , w„ + 1) and v = c "~"'^ * (wi — n, . . . , w„) . At most 
two monomials get added to M (v is added only if its coefficient is nonzero). Wc 
remove the first monomial (suppose this is u) and apply the reduction again. So 
U2 = cA^ * {wi - n + 2,W2 - n + 2, . . . ,Wn + 2) and % = cA ""^"'^ * {wi - n + 
1, W2 — n + 1, . . . ,Wn + 1). Once more, at most two monomials get added to M. 
But now notice that the next item on the list would be v. Applying the reduction 
to this monomial would give W3 = cA ""^"'^ * {wi — n + 1,W2 — n + 1, . . . ,Wn + I) 
and V3 = c ""^™^ * — n,W2 — n, . . . , w„). Notice that W3 = kv2, where k is 
some constant, and thus for each pair of monomials we remove we are only really 
adding, at most, the number of monomials we had before plus one. 
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One can visualize the process in a diagram as follows: 

{wi,W2, . . . , Wn) 




Wn + 1) 



(n — wi ) /n 

(wi — n, W2, . 

(n — W2)/n 

{wi — n, W2 — n, . . . , Wn) (wi — n + 1,W2 — n + 1, . ■ ■ , w„) {wi — n + 2, W2 — n + 2, . . . , Wn + 2) 

Notice the down arrows indicate subtracting n from the first position that is 
greater than n, and this process clearly terminates. Also notice that the "down 
and right" diagonals involve subtracting n from the first entry larger than n and 
then adding 1 to all entries. If this process does not terminate in a monomial whose 
entries are all less than n then at the n-th reduction step we will find the monomial 
{wi — n + W2 — n + n, . . . , Wn — n + n) — w, and as long as we can empty the list 
M , the algorithm terminates. But since all of the monomials are accounted for by 
following the "down and right" diagonal and the down arrows, we are certain we 
can empty M . □ 

The following proposition is a direct consequence of Algorithm [T] 

Proposition 3.1. W is generated over C(A) hy the set of monomials 

B = ==a;'"|l < w, <n-l,^uH = mod n| . 

In particular, W has dimension 

[n - ~{n~ + (n _ i)(«-3) ± (n - 1). 

The vector bundle W is equipped with an integrable connection V defined by 



d_ 
dX 



In particular, for monic monomials, 

OA 

Katz, in |16', proved the following useful lemmas. 

Theorem 3.2 (The Comparison Theorem). Let wq = ■^J27=i''^i- There is a T- 
linear map 7^ : ^ iJ^^'^^P" - Xx/T) given by 



By this theorem and the residue map (cf. [13 ) we have an isomorphism from 



H 



dR 



-(P«-Xa/T) toH^-'{Xx/T) 



And we have the following: 
Theorem 3.3. The map Q induces, by passage to quotients, an isomorphism 

Q:W^H2r\XxIT), 
which is compatible with the connection. 
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Hence the space W obtained through Dwork's construction is isomorphic to tlie 
middle (relative) cleRham cohomology. 

It is also important to note that Q transforms V into the Gauss-Manin connec- 
tion. 

4. Computing the connection matrix 

Let /LtJ^ denote the group of n-tuples of n-th roots of unity, and A denote the 
diagonal elements. 

The character group of /U^/A is in one-to-one correspondence with the set 
W = {{wi, . . . ,Wn)\0 <Wi<n, ^ = mod n}, 

where 

x«,(0:=x(r), r = c---c 

and X is a fixed primitive character of 
Let 

G = {ee<|a---?n = i}/A. 
The characters Xw of A*" /A which act trivially on G are precisely powers of xy> 
where 1 = (1,1,..., 1). Thus, Char(G'), the character group of G, corresponds to 
equivalence classes of w in W, where w' w \i w — w' is a, multiple (mod n) of 1. 

The varieties Xx allow a faithful action of the group G by ^ = (^i , . . . , taking 
the point {xi,.. . ,a;„) to (^ia;i, . . . ,£,nXn)- Using this action, we get that W splits 
into eigenspaces as follows. 

Proposition 4.1. The action of G on a fiber W, gives 

w = W(x), 

XGChar(G) 

where W(x) is an eigenspace with basis 



{w, w + l(modn), . . . ,w + n — l(modr?.)}. 
but we exclude adding any vector m such that m + Wi = mod n for some i. 

To understand V's effect on W, it is sufficient to know what it does to elements 
in the basis B. Prom the definition of V we see that 

V(x«') = ■^Fxx'" = -nx"'+T 

where w + 1 ~ (wi + 1, . . . , w„ + 1). 

Thus, the connection commutes with the action of G, so the proposition implies 
that V preserves eigenspaces. We want to compute the connection matrix V. 
Because of the way in which V preserves eigenspaces, the connection matrix will 
have blocks on its diagonal for each eigenspace. 

The main idea of the following algorithm is to use the reduction algorithm de- 
scribed earlier on V(a;"') where w runs through the basis of an eigenspace. 

Algorithm 2 (Computing a block of the connection matrix). This algorithm takes 
any vector of integers as an input and outputs a matrix that is the block of the 
connection matrix that corresponds to that vector's eigenspace generators. 
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(1) Create a basis of the eigenspace related to w by computing B = {w,w + 
l(modn), . . . , w + n — l(niodn)} = {wi, . . . , Vk}, where we omit monomi- 
als which have entries equal to mod n (so k may or may not equal n). 

(2) Create M , a k x k matrix of zeros. 

(3) Let i = 1. 

(4) If i = k-\-\, output M . Otherwise, take the monomial Vi in B and compute 
its derivative, that is V(ui) = —n * {vi + 1). 

(5) // all the entries in V(t;i) are less than n, then add —n to the {i,i + 1) 
position in M, add 1 to i, and go back to step 5. If not, proceed to the next 
step. 

(6) Reduce V(wi) according to Algorithm{l\ From the algorithm, we have V(wi) 
as a linear combination of elements of B, so that V(wi) = aiVi + ■ ■ ■ akVk. 
Replace the i-th row of M by (ai, . . . , ak). Add 1 to i and go to Step 5. 

5. The differential equation associated to the connection 

In this section, we will show that the differential equation associated to the 
connection V is a hypergeometric differential equation. We have developed an 
algorithm which outputs the parameters a, /3 given each block representative. 

We will first establish some notation and definitions. For more details and proofs, 
see for example [3], [5], or [T^ . 

Consider the nth order equation 

n 

(1) an-m(^)2/^'"^ = 0, (ao(z) = 1) 

m=0 

where the ak{z) are single- valued and analytic in a punctured neighborhood of a 
point zq. Recall that if any of the ak have a singularity at zq, then zq is called a 
singular point for ([T]), otherwise it is called an analytic point. We say zq is a regular 
singular point if 

(2) ak{z) ^ {z - zoy''bk{z), (/c = 1, . . . , n), 

where b^ is analytic at zq. 

A system of n first order equations over C(z) has the form 

(3) y' = Ay 

in the unknown column vector y = (j/i, . . . , 2/„)'^ and where A is an n x n-matrix 
with entries in C{z). The entries are assumed to be single-valued and analytic at a 
neighborhood of a point zq, and will at most have a pole at that point. 

If A has a singularity at zq, then zq is a singular point for the system ([3]). zq is 
a regular singular point if 

A{z)^{z^zo)-^A{z) 

where A is analytic for a neighborhood of zq (including zq), and A{zq) ^ 0. 

To study the system at oo, we change variables from z to l/C- The associated 
system is 
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It is not difficult to see that a differential equation like ([T]) can be rewritten 
as a system by setting yi = y,y2 = y' , ■ ■ ■ ,yn = y^""^-*. Notice that this means 
y'l = 2/2,2/2 = 2/3: ■• ■ ,2/n-i = 2/rn ^-^d y'„ is given by the differential equation. So 
the differential system is determined by a companion matrix, as follows: 





1 \ 




( 





1 









2/1 \ 


d 


2/2 












1 







2/2 


dz 






















\ 2/n / 






-a„ 




-a»i-2 ■ • 


■ -ai / 


V 


2/n / 



(4) 



Now, we will explain how a system of first order differential equations arises from 
the connection. Recall that, on a vector bundle, being equipped with a connection 

V is equivalent to being equipped with an action of — (see for example pT] for a 

d\ 

detailed explanation). In short, we have a first-order system defined by 



dxy = ^y\ 

where A is actually the transpose of the matrix we found in the previous section. 
In fact, each block defines its own differential system. We will use the following 
fact: 

Theorem 5.1 (Cyclic Vector Lemma). Any system oj linear first order differential 
equations is equivalent to a system which comes from a differential equation. 

A proof of the lemma can be found in [T^] . 

Basically, this theorem says that in the space of solutions of our differential 
system there is a cyclic vector, that is, a vector such that v, Av, A'^v, . . . , A"^^v 
is a basis. More specifically, this means that if y = (yi, 7/2, • • ■ , 2/n)^ is a solution 
for the system, we can find an equivalent system with solutions of the form y ~ 
dyi d"-!"- 

j/„ . This determines a change of basis matrix 
S such that Sy = y. 

The vector Sy = y satisfies a differential system of the form 



2/1, 
(fc) 



) . In fact, using the system, we can represent the derivatives 
y-j^ ' as a linear combination of yi , 



y, 



and this last system is the companion matrix to a higher order differential equation. 

In our situation, since the basis vectors are basically already derivatives of each 

other, any vector in the basis, for example yi, is a cyclic vector, and so S is easy 

/ 

to determine. Let As = ( SAS^^ + —S' 

\ d\ 

Notice that, from the reduction algorithm, the entries in the connection matrix 
are polynomials in A or rational functions in A. In the case that they are rational 
functions, 1 — A" is the only possible denominator. After doing the change of basis 
described above, we may get some powers of A in the denominator, but as the new 
system is a companion matrix, this will only happen in the last row of the matrix. 

At this point, we have a differential equation associated to the connection. But 
solving high order differential equations is not a simple task. Instead, we will show 
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that these matrices are related to a hypergeometric group, which in turn gives us 
the defining parameters of the hypergeometric differential equation we want to find. 

We have a way of changing from a differential equation to a system and viceversa. 
It is important to note that a regular singular point of ([ij zq may not be a regular 
singular point of the system associated with it. This happens only when the Uk 
have at most simple poles at zq. However, there is an equivalent first-order system 
with the property that if zq is a regular singular point of ([T]) then zq is a regular 
singular point of the system. 

Suppose H]) has a regular singularity at zq, and let be a solution of ([1]). Define 
(j) to be the vector with components 0i , . . . , 0„ by setting 



Then, clearly, 



(z-zo)^', = iz-zo){{z-zo)'-'cj,<^'-'^y 

= {z- zo){{k - l)(z - zo)'^"'0('='') + iz~ zo)'-V^'^) 
= (k- l)(z - zo)"-'^^'-''' + iz- Zo)'(b^''^ 
= (fc - l)0fc + 0fc+i (fc = l...,n-l) 

And, finally, 

n 

{z - Zo)(/)'„ = {n- l)(l>n - ^ bn-,n+l{z)(f>m, 

rn—1 

where the bi are defined as in ([2]). 

Therefore (p is a. solution of the linear system 

(5) y' = A{z)y 

where A has the structure 



A{z) = (z-zoY 



( 







V -6„ 



1 
1 






-bn- 




1 
2 




-K-2 





1 
3 





(n-l)-fei / 



This system clearly has a regular singularity at zq. 

In our situation, we want a hypergeometric differential system, so in particular 
we want a simple pole at A = 0. We can accomplish this by replacing the matrix 
As by the matrix A, which has a simple pole at 0. In particular, now we are certain 
that the only denominators in the last row are of the form 1 — A" . 

We will now change variables to z = A" . By the chain rule, we have: 



d d d d 

-ttV = -rV^^ TV * "^^ 
d\ dz d\ dz 



-1) 
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So our system representing the derivative with respect to z is y' — By where 

B ~ — — T^A. Since A has a simple pole at A = 0, replacing every instance of A" 

by z gives that B has a simple pole at 2 = as well, but now we also have a simple 
pole at z = 1. 

There is an algorithm by Brieskorn which relates Gauss-Manin connections to 
monodromy group generators [3]- Let A be the matrix representation of the con- 
nection. The algorithm uses the fact that if A has a simple pole around a given 
point, i.e. can be written as 

A = A_i(z - zo)-^ + Ao + Ai(z - zo) + • • • 

then R = e'^'^^^-^ gives the monodromy around zq. 

The process of changing the system matrix to B ensures that we have a simple 
pole around zero, one, and infinity. Since this last system has regular singular 
points at 0, 1, 00 and no other singularities, it is Fuchsian, as we expected. 

Let B-i denote the residue around zero, and B^i denote the residue at 00. Now 
define ho ~ g^^TviB^i ^^^^ _ g27rjB_i j^g^- jj residue of B around z ~ 1. 

The matrix hi = e^'^*^ is clearly a reflection in the sense described by Beukers 
and Heckman. This is easily checked by noticing that there is only one row with 
denominators of the form 1 — z. Thus, the residue at one will necessarily have rank 
one, which implies that hi — Id has rank one. 

We have shown that the matrices hoc, hi, /iq generate a hypergeometric group. 

Recall that the monodromy group described by Beukers and Heckman in [2] 
is generated by the monodromy matrices around zero, one, and infinity, and the 
parameters are determined by the eigenvalues of these matrices. In fact, we can 
get the parameters of the hypergeometric differential equation directly from i?_i, 
the residue around zero, and -B-i, where B is the system at cxo. The eigenvalues 
of B-i will be the /3's and the eigenvalues of B-i will be the a's. Therefore, the 
monodromy group we have found corresponds to the differential equation 



D{a; f3\z)y = 0, 

where the a's are the eigenvalues of and the /3's are the eigenvalues of B-i. 
Here is the algorithm we have just described. 

Algorithm 3 (Computing the parameters of the hypergeometric differential equa- 
tion). This algorithm takes as input a monomial w and computes the parameters 
of the hypergeometric differential equation associated to the eigenspace generated by 
w. The output is two lists, a and P, containing the parameters. 

(1) Let A be the k x k connection matrix block associated to w, computed with 
Algorithm\^ 

(2) Compute the change of basis matrix, S , given by assuming w as a cyclic 
vector. 

(3) Compute As = (^SAS'^ + ^S-^^ . 

(4) Compute A by multiplying Ukj by A''^-'^^ and adding i — 1 in the {i,i) 
position of As for i — I, . . . ,k . This means that in the (fc, k) position we 
have k ~ 1 + Xak.k- 

(5) Multiply A by l/n, call this B. 
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(6) Let ho = Resz=oB and compute the eigenvalues. Let /3 be the list of eigen- 
values. 

(7) Let hoo = ReSz=ooL3 and compute the eigenvalues. Let a he the list of these 
eigenvalues. 

(8) Output: a, /3. 

Remark 5.1. While computing some examples, we noticed that given a vector 
(wi, W2, . ■ ■ , Wn), if we cancel out the numbers which it has in common with the 
list (0,1,2,. — 1), then at — ^ for each Wj that survives the cancelation, and 
Pi — ^ for each k that survives in the second vector. This is exactly Katz's main 
result in the case of hypergeometric sheaves [18] . 



6. An Illustrative Example 
Suppose n = 6. By Algorithm [T] we have that 



W = (x""]! <w, <5,^u;, =0mod6^, 

and W has dimension 5^ - S'' + 5^ - 5^ + 5 = 2605. 

Take w — (1, 1, 1, 2, 2, 5). This belongs to the eigenspace generated by 



%,i,i,2,2,5) = {(1,1,1,2,2,5),(3,3,3,4,4,1),(4,4,4,5,5,2)}. 

Here is an example of the algorithm for computing the block in the matrix 
representation of V for n = 6 and the eigenspace corresponding to the monomial 
(1,1,1,2,2,5) ~ xiX2X^x\x\x\, with basis denoted earlier by i i_2,2,5)- I will 
denote this block by Vb^i 11225)- 

(1) Apply V(l,l,l,2,2,5) = -6(2,2,2,3,3,6). Using the relations we can 
write this last monomial in terms of the monomials in ^(1^1^1^2,2.5)- We can 
represent the process of Algorithm [1] graphically, as shown below: 



(2,2,2,3,3,6) 




(2,2,2,3,3,0) (3,3,3,4,4,1) 

This means that (2, 2, 2, 3, 3, 6) = A(3, 3, 3, 4, 4, 1) + - (2, 2, 2, 3, 3, 0), which 
is a monomial in ;S(i^i_i^2.2,5)- 

Thus, in the matrix representation of Vg^^ 11225)' there will be a — 6A 
as the (2,1) entry. 

(2) We repeat this process for the next monomial in the basis, (3, 3, 3, 4, 4, 1). 
Applying the connection we get V(3, 3, 3, 4, 4, 1) = —6(4, 4, 4, 5, 5, 2). Since 
this monomial is already in B we write —6 in the (3, 2) position in the block 
matrix. 

(3) Take V(4, 4, 4, 5, 5, 2) = -6(5,5,5,6,6,3). We have to do the reduction 
process again, represented below. 
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(5,5,5,6,6,3) 



(5,5,5,0,6,3) (6, 
-1/6 



,1,7, 4) 



(6,6,6,1,1,4) 



(7, 7, 7, 2, 2, 5} 
-1/6 / \ 



(0,6,6, 1, 1,4) 



(1,7,7,2,2,5) (2,8,8,3,3, 
-1/6,/ ^\ \ -2/6 / \ X 



(1,1,7,2,2,5) (2,2,8,3,3,6) (3,3,9,4,4,7) 

-1/6 X ^\ -2/6 ^\ -3/6 



(1, 1, 1, 2, 2, 5) 



(2,2,2,3,3,6) (3,3,3,4,4,7) 
/ \ ^ -1/6 , 



(4,4,4,5,5, 



-2/6 



(2,2,2,3,3,0) (3,3,3,4,4,1) (4,4,4,5,5,2) (5,5,5,6,6,3) 

This is a bit harder to unravel than the other cases, but it works in 
exactly the same way. The diagram shows us that 



\2 i'T\4 

(5,5,5,6,6,3) = - — (1,1,1,2,2,5)+ ^(3,3,3,4,4,1) 
3 \ 5 

- — (4,4, 4,5,5,2) + A'^(5,5,5, 6,6,3). 
And solving for (5, 5, 5, 6, 6, 3) we get that 

V(4, 4, 4, 5, 5, 2) = -6(5, 5, 5, 6, 6, 3) 

= -T8(A^(^' ^' ^' ^' + W^)^'' 



9A5 



A6-1 



(4,4,4,5,5,2). 



(4) Combining all of these steps, we can write Vb^^ ^ i ^ s ^, as 



^B(i,i.i 



(1,1,1,2,2,5) 



V 




-6A 

-6 



A2 



18(A6 - 1) 

17A^ 
6(A6 - 1) 
9A^ 

~A6^ / 



We will now describe the algorithm for finding the parameters of the differential 
equation for the same example monomial. The connection above gives us that 
differentiation with respect to A is equivalent to the differential system 
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-6A 





\ 


d 










-6 










17 A* 


9A5 






V 


18(A6 - 1) 


6(A6 - 1) 


A6-1 


/ 



We compute the change of basis for the cychc vector lemma, which is 



S = 



10 
-6A 
-6 36A 



And so we have the new system 



4.y=isAS-^ + '4s 



dX 



dX 



y 






2A3 



lOA^ 



7A6 



y, 



V 1-A6 A2(l-A6) A(l-A6) / 

which, as we expected, is given by the companion matrix for a third order differential 
equation. 

This system is equivalent to 



_ 1 
dx'^~ X 



( ° 



2A6 
V T^A6 



1 
1 

lOA^ 





1 

7A6 



y, 



A6 



1-A6 

which clearly has a simple pole at A = 0. 

Now we can change variables by setting z = X^. The change of variables leaves 
us with a system 



dX z 



The residue at zero is 



1/6 
1/6 
bz-l 





1/6 
5^ + 4 



\ 3(1 -^) 3(1-0) 6(1 -^) 



which has eigenvalues 0, 1/2, 1/3. 
Also, around infinity we have 



1/6 
1/6 1/6 
-1/3 2/3 



dC 



-1/6 
-1/6 

5-C 





-1/6 
5 + 4C 



V 3(1-0 3(1-0 6(1-0 / 
which has residue (at ^ = 0) of 

-1/6 

-1/6 -1/6 

1/3 5/3 5/6 
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and thus yields the eigenvalues 1/3, 1/6, 1/6. 

We have now found the parameters of the hypergeometric differential equation 
associated to this connection matrix block: 

^ /I 1 1 1 2\ 

To sum it up, the block of the matrix V corresponding to the eigenspace of 
(1, 1, 1, 2, 2, 5) gives rise to the hypergeometric differential equation which has 



3F2 



11112 



6' 6' 3' 2' 3 
as its holomorphic solution around 0. 

Table 1 below shows some numerical examples for n = 6. 



Table 1. Parameters for n = 6 



Monomial 




ft 


[1,1,1,1,1,1] 


[11111 

6 ' 6 ' 6 ' 6 ' 6 






12 5 1 

2' 3' 6' 3 




[5,3,1,1,1,1] 




"1 1 1 

6 ' 6 ' 6 






2 1 
3' 3 




[4,4,1,1,1,1] 




1112" 

6' 6' 6' 3 






"1 1 5" 

3' 2' 6 




[5.'?. 2. 1.1.1] 




"1 1 1 
3' 6' (j 






1 2 
2' 3 




[4,3,2,1,1,1] 




1 1" 

6' 6 






"5" 

6 




[3,3,3,1,1,1] 




1111" 
2' 2' 6' 6 






"2 5 1" 
3' 6' 3 




[4,2,2,2,1,1] 




"1 1 1 

3' 3'6 






1 5 
2' 6 




[3,3,2,2,1,1] 




"111 
2' 3' 6 






2 5 
3' 6 




[3,2,2,2,2,1] 




"1 1 1 

3' 3'3 






2 5 
3' 6 




[5,5,3,3,1,1] 




"1 1 5 

6' 2'6 






1 2 

3' 3 




[5,5,4,2,1,1] 




1 5" 

6' 6 






"1" 

2 




[5,4,4,3,1,1] 




1 2" 
6' 3 






"1" 
3 




[5,4,3,3,2,1] 




"11 

2 


D 


[4,4,4,3,2,1] 




2 2" 
3' 3 






"5" 
6 
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Appendix 
GP Scripts 

Computing the connection matrix. This function counts the number of coor- 
dinates of a vector with entries bigger than or equal to d. 

count (a) = 
{ 

local(t) ; 
t=0; 

for(k=l,length(a) , if (a[k] >=length(a) ,t=t+l)) ; 
t 

} 

This function (from [22j) can tell if a given element is in a vector, and gives the 
"position" of the element. 

memb(g,v)=f or (k=l , length (v) ,if (g==v[k] ,return(k))) ;0 

The following function takes a vector a and an integer m (its coefficient) and 
subtracts one to all the entries and adds d to one of them until it gets to 0, a vector 
in the basis, or the original vector. It saves the leftovers in a vector v. It is one of 
the two possible reductions coming from the relations on W. 

redl (a,m)= 
{ 

local (j ,l,b,h,t,d,s,v,u) ; 

h=m; 

b=0; 

t=vector (length(a) ) ; 

d=length(a) ; 

1=1; 

j=i; 

s=vector (length(a) ) ; 

v=vector (0) ; 

t=a; 

if (count (t)>0, until (count (t)==0 II t==a, 

for(k=l,length(t) , if (t [k] <t [j] , j=k)); 
for(k=l,length(t) ,s[k]=t [k]-l) ; 

for(k=l,length(t) , if (k==j , t [k] =s [k] +d, t[k]=s[k])); 
l=h*s [j] /(d*n) ; 
v=concat (v , [ [1 , s] ] ) ; 
h=h*l/n)) ; 
[h,t] 

} 



redllef tovers (a,m)= 

{ locaKj ,l,b,h,t,d,s,v,u) ; 

h=m; 

b=0; 

t=vector (length(a) ) ; 
d=length(a) ; 
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1=1; 

j=i; 

s=vector(length(a)) ; 

v=vector (0) ; 

t=a; 

if (count (t)>0, until (count (t)==0 II t==a, 

for(k=l,length(t) , if (t [k] <t [j] , j=k)); 
for(k=l,lengtli(t) ,s[k] =t [k]-l) ; 

for(k=l,length(t), if (k==j , t [k] =s [k] +d, t[k]=s[k])); 
l=h*s[j]/(d*n) ; 
v=concat ( v , [ [1 , s] ] ) ; 
h=h*l/n)) ; 

V 

} 

Here is the other possible reduction. This one subtracts n from one spot and 
adds one to everything afterwards. Saves leftovers in vector v. 

red2(a,m)= 
{ 

local (j ,l,b,h,t,d,s,v,u); 

h=m; 

b=0; 

d=length(a) ; 
1=1; 

t=vector(length(a)) ; 
v=vector(0) ; 

j=i; 

s=vector (length (a)) ; 
t=a; 

if (count(t)>0,until(coiuit(t)==0 II t==a, 

for(k=l,length(t) , if (t [j] <t [k] , j=k)); 

for(k=l,length(t) ,if (k==j ,s[k]=t [k] -d,s[k]=t [k] ) ) ; 

f or (k=l, length (t) , t[k]=s[k]+l); 

l=h*(-s[j])/d; 

v=concat (v , [ [1 , s] ] ) ; 

h=h*n) ) ; 

[h,t] 

} 

red21ef tovers(a,m)= 
{ 

local (j ,l,b,h,t,d,s,v,u); 

h=iii; 

b=0; 
d=length(a) ; 
1=1; 

t=vector (length(a) ) ; 
v=vector (0) ; 

j=i; 
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s=vector (length (a) ) ; 
t=a; 

if (count(t)>0,until(count(t)==0 II t==a, 

for(k=l,length(t) , if (t [j] <t [k] , j=k)); 

for (k=l , length (t ) , if (k== j , s [k] =t [k] -d, s [k] =t [k] ) ) ; 

for(k=l,length(t) , t [k] =s [k] +1) ; 

l=h*(-s[j])/d; 

v=concat (v , [ [1 , s] ] ) ; 

h=h*n) ) ; 

V 

} 

Now we combine these two reductions and loop until we get monomials in the 
basis of W. The input of this function is a vector of any length and the output will 
be the "linear combination" of that vector in terms of the basis vectors (vectors 
with entries between 1 and the length). 

reduction(a)= 
{ 

local(d,b,c,u,v,w,uu, bb,j, t, s, g,r); 
u=vector (0) ; 
v=vector (0) ; 
d=length(a) ; 

j=i; 

if (count (a)==d, u=[redl(a,l)] ; v=redllef tovers (a, 1) , 

u=[red2(a,l)] ; v=red21ef tovers (a, 1)) ; 

for(k=l, lO'd, 

if (k<=length(v) , 

if (count (v[k] [2])==0, 

if (v[k] [1]==0, ,b=0; 
for(i=l,length(u) , 

if(v[k] [2]==u[i] [2], 

u[i] [l]=u[i] [l]+v[k] [1], 
b=b+l)); 

if (b==length(u) ,u=concat (u, [v [k] ] ) ) ) , 
if (v[k] [1]==0, , 

uu=red2(v[k] [2] ,v[k] [1] ) ; 

v=concat(v,red21ef tovers (v[k] [2] ,v[k] [1])) ; 
b=0; 

for(i=l,length(u) , 

if (uu[2]==u[i] [2] , 

u[i] [l]=u[i] [l]+uu[l] , 
b=b+l)); 

if (b==length(u) ,u=concat (u, [uu] ) ) ) ) , 

break) ) ; 

b=0; 

for(k=l, length(u), 
if (u[k] [2]==a, 

w=vector (length(u)-l) ; 

for(j=l,k-l,w[j] = [u[j] [l]/(l-u[k] [l]),u[j] [2]]); 
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f or ( j =k , length (u) - 1 , 

w[j] = [u[j+l] [l]/(l-u[k] [l]),u[j+l] [2]]), 
b=b+l)); 

if (b==length(u) , r=u, r=w) ; 

r 

} 

The next step is to write the connection matrix from this, that is, write a function 
that gives the derivatives of each vector in terms of the basis. In fact, there is an 
easy way to write the derivative of any vector using the reduction function. 

derivative (a) = 
{ 

local(d,t,w) ; 

d=length(a) ; 
t=vector (d) ; 
for(k=l,d,t[k]=a[k]+l) ; 
if (count (t)==0, [[-d,t]] , 

w=reduction(t) ; 

for(k=l,length(w) , 

w[k] [l]=w[k] [l]*(-d)); w) 

} 

Given a basis vector, we can find all the other basis vectors that will be a basis 
for the same eigenspace. 

orbit (a) = 

{ local(l,iii, c,ss); 
d=length(a) ; 
1=0; 
c=0; 
ss=l; 

for(k=l,d, for(t=l,d, if (k==d-a[t] , 1=1+1 ; break) )) ; 
ni=d-l ; 

b=vector(in) ; for(k=l,ni, b [k] =vector (d) ) ; 
b[l]=a; 

f or(s=l,d-l,f or(t=l,d, if (s==d-a[t] , , c=c+l)); 

if (c==d,ss=ss+l;for(t=l,d, b[ss] [t] = (a[t]+s)7.d)) ;c=0) ; 

b; 

} 

The following gives the matrix representation of the block of the Gauss-Manin 
connection associated to a particular basis vector (i.e., it gives a block of the whole 
matrix, which is related to the eigenspace related to this basis vector). 

connectionmatrix(a)= 
{ local (v,w,M) ; 
v=orbit (a) ; 

M=matrix (length (v) , length (v)) ; 
f or(j=l,length(v) , 

w=derivative (v [ j ] ) ; 

f or (k=l, length (w) ,M[memb(w[k] [2] ,v) ,j]=w[k] [1])) ; 
M=mattranspose(M) ; 
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M 

} 

The algorithm to find the differential equation. The following finds the de- 
rivative with respect to A of a vector with a coefficient. It is essentially the product 
rule. 

derivn(a)= 

{ local (b,z, WW, vv) ; 
b=deriv(a[l] ) ; 
z=derivative (a [2] ) ; 
ww=[[b,a[2]]] ; 

vv=vector (length(z) ) ; 

f or (k=l, length (z), vv[k] = [a[l]*z[k] [1] ,z[k] [2]]); 
for(k=l,length(vv) , 

if (vv[k] [2]==a[2] , ww[l] [l]=ww[l] [l]+vv[k] [1] , 
ww=concat (ww , [vv [k] ] ) ) ) ; 

WW 

} 

We would like to have the derivative of a vector which is a linear combination of 
these monomials. This should use ideas like the function above. The first function 
finds the derivative of a vector (with a coefficient) in a prescribed basis determined 
by the orbit of b. The second does the same, but only outputs the vector of 
coordinates, without writing the basis down. 

derivv(a,b)= 
{ local (v,w,z) ; 
w=orbit (b) ; 
v=vector (length (w)) ; 
for(i=l,length(w) , v[i]=[0,w[i]] ) ; 
for(k=l,length(a) , 
z=derivn(a[k] ) ; 
for(j=l,length(z) , 

v[memb(z[j] [2] ,w)] [l]=v[memb(z[j] [2] ,w)] [l]+z[j] [1] 

); 

); 

V 

> 

derivv2(a,b)= 
{ local (v, w, z) ; 
w=orbit (b) ; 
v=vector (length (w)) ; 
for(k=l,length(a) , 
z=derivn(a[k] ) ; 
f or(j=l,length(z) , 

v[memb(z[j] [2] ,w)]=v[meinb(z[j] [2] ,w)]+z[j] [1] 

); 

); 

V 
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} 

We want to change basis, and we need a matrix that changes from our basis 
obtained by using "connection" to a basis obtained from derivatives (i.e. we use the 
cychc vector theorem to write the matrix as the companion matrix to a differential 
equation). 

cob(a)= 
{ 

local (z, vv, uu, w) ; 
uu=orbit (a) ; 

r=vector (length (uu) ,k, if (k==l , 1) ) ; 
vv=vector (length (uu) , k, [r [k] ,uu[k]] ) ; 
z=[r] ; 

for (k=l , length (uu) -1 , 

w=deri vv2 ( vv , a) ; 

vv=derivv(vv,a) ; 

z=concat (z , [w] ) 

); 

S=Mat(z~) ; 
S 

> 

This function takes two input vectors, one is a basis and the other a vector 
indicating a hnear combination of elements in this basis. 

cobv(uu,r)= 
{ 

local (z, vv, w) ; 

vv=vector (length (uu) , k, [r [k] ,uu[k]] ) ; 
z=[r] ; 

for (k=l , length (uu) -1 , 

w=derivv2 (vv , a) ; 

vv=derivv(vv,a) ; 

z=concat (z , [w] ) 

); 

S=Mat(z~) ; 
S 

} 

The following computes what a change of basis does to the system of differential 
equations, where we change from a basis found by using the connection function to 
a basis of all the derivatives of a specific vector. 

cobsystem(A, S)= 
{ 

local (dS, C) ; 

dS=matrix (length (A) , length (A) , X, Y, deriv(S [X, Y] ) ) ; 
C=S*A*l/S+dS*l/S; 
C 

} 

Now, as seen in Section [SJ we can change this system into an equivalent one with 
a simple pole at 0. 
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regf orm(A)= 

{ 

local (m.Areg) ; 
in=length(A) ; 
Areg=A; 

f or (k=l ,m, Areg[m,k] =Areg[m,k] *n" (m-k+1) ) ; 
f or (k=l ,in, 

for(i=l,m, if (i==k, Areg [k, i] =Areg [k, i] +k-l) ) ) ; 
Areg 
} 

We need to change variables to have it in terms of z instead of A**. 

varchange ( A , d) = 
{ 

for(k=l,length(A) , 

for(i=l,length(A) , 

A[k,i]=substpol(A[k,i] ,n"d,x)/d)) ; 

A 

} 

Finally we can compute the residue of the matrix at z = ^. 
residuezero(A)= 
{ 

for(k=l,length(A) , 
for(i=l,length(A) , 
A[k,i]=subst(A[k,i] ,x,0))) ; 
A 
} 

We want a function that finds the rational roots of a polynomial with rational 
coefficients (because all of our characteristic polynomials are of that form and only 
have rational roots). We are using the rational roots theorem. 

ratlroots(f )= 
{ 

local(p,q,z,vv,a,b,c,n,r, j ,i) ; 

c=poldegree(f ) ; 
vv=vector (c+1) ; 
r=vector(0) ; 
z=vector (0) ; 

for(k=l,c+l, vv[k]=polcoeff (f ,k-l)) ; 
n=denominator (vv) ; 
f =f *n; 

a=substpol(f ,x,0) ; 

if(a==0, r=concat(r,0) ; f=f/x; a=substpol(f ,x,0)) ; 

b=pollead(f ) ; 

p=concat (divisors (a) , -divisors (a) ) ; 
q=concat (divisors (b) , -divisors (b) ) ; 
f or (k=l, length (p) , 

f or(i=l,length(q) , 

if (memb(p[k]/q[i] ,z)==0, z=concat (z,p [k] /q[i] ) ) ) ) ; 
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for(i=l,length(z) , 

if (substpol (f ,x, z [i] )==0, r=concat (r ,z [i] ) ; 

f=f/(x-z[i]);i=l)); 

r 

> 

We now want to combine all these steps to find the hypergeometric parameters 
given a vector (or monomial) in W. 

hypergcoef f (a)= 
{ 

local(A, B, S, Azero, Ainf,f,g, vv, uu, d, m,r,t); 

d=length(a) ; 

A=connectionmatrix(a) ; 

S=cob(a) ; 

A=cobsystem(A , S) ; 

A=regf orm(A) ; 

A=varchctnge (A , d) ; 

Azero=residuezero (A) ; 

f=charpoly (Azero) ; 
B=matrix (length (A) .length (A)) ; 
Ainf =matrix(length(A) ,length(A)) ; 
f or(k=l ,length(A) , 

f or(i=l ,length(A) , 

B[k,i]=-substpol(A[k,i] ,x,l/y))) ; 
f or(k=l ,length(A) , 

f or(i=l ,length(A) , 

Ainf [k,i]=subst(B[k,i] ,y,0))) ; 

g=charpoly(Ainf ) ; 
r=ratlroots(f ) ; 
m=vector (length (r) ) ; 
for(k=l,length(r) , m [k] =l-r [k] ) ; 
\\ printC'alphas " ,ratlroots(g) , " betas ",m) 
t= [ratlroots (g) ,m] ; 
t 

} 

Creating a hypergeometric table. To make a table like Table IH] we should be 
able to check the hypergeometric coefficients somewhat systematically. So we have 
to find a good way to generate basis vectors (or representatives up to permutations 
of the variables). We are certain that there are more efficient ways to do this, based 
on conversations with computational number theorists, but the focus of our project 
is to obtain the relationship with hypergeometric functions, which this code does. 
This function (also from [22] ) finds all the partitions of a number m. 

part(m)= 
{ 

local (k, j ,sm,sj ,s, S = [] ) ; 
k = j = 1; 

sm = sj = vector (m+1) ; 
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while (k, 

s = sm[k] +j ; 
if (s > m, 

untilCj <= m, j = sj[k]+l; k — ); 
next) ; 

k++; sni[k]=s; sj[k]=j; 
if (s == m, 

S = concat(S, [vector (k-1,1, sj [k-1+1] )] ) ) ) ; 

S 

} 

This function uses the previous one to find the partitions of a number m of length 
c and into numbers that are less than c. 

part2(m, c)= 
{ 

local (v,w,t) ; 
v=part (m) ; 
w=vector (0) ; 
f or(k=l ,length(v) , 

if (length (v[k] )==c, 
t=count (v [k] ) ; 
if (t==0, 

w=concat (w , [v [k] ])))); 

w 

} 

Now we can put this together to get representatives of the basis. We don't 
have strict representatives, but at least we eliminate the cases in which none of the 
entries are equal to one, because those obviously are in an eigenspace with a vector 
with entries equal to one. 

basisreps(m)= 
{ 

local (v,w) ; 
v=vector (0) ; 
for(k=l,ceil((m-l)/2) , 

w=part2(m*k,m) ; 

for(j=l, length (w) , 

if (memb ( 1 , w [ j ] ) ==0 , , v=concat (v,[w[j ]])))); 

V 

} 

We now want to be able to output a table with the basis vectors and the hyper- 
geometric parameters associated to them given a number d. It should turn out to 
be the list of numbers in the vector that remain after canceling out with the list of 
numbers between and d, in accordance with Katz's results [TS] . 

hypergtable (d) = 
{ 
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local (v,u) ; 
v=basisreps (d) ; 

f or(k=l ,length(v) , u=hypergcoef f (v[k] ) ; 
print (v[k]," ", "alphas ", u[l], 

betas ", u[2])) 

> 
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